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We study the thermal conductivity of the disordered two-dimensional electron gas. To this end we 
analyze the heat density-heat density correlation function concentrating on the scattering processes 
induced by the Coulomb interaction in the sub-temperature energy range. These scattering processes 
are at the origin of logarithmic corrections violating the Wiedemann-Franz law. Special care is 
devoted to the definition of the heat density in the presence of the long-range Coulomb interaction. 

To clarify the structure of the correlation function, we present details of a perturbative calculation. 
While the conservation of energy strongly constrains the general form of the heat density-heat density 
correlation function, the balance of various terms turns out to be rather different from that for the 


correlation functions of other conserved quantities 
density correlation function. 

PACS numbers: 71.10.Ay, 72.10.-d, 72.15.Eb, 73.23.-b 

I. INTRODUCTION 

While experimentally thermal transport is controlled 
by boundary conditions, for the theoretical description it 
is more convenient to study the response to a gradient 
of temperature. A principle difficulty for the description 
of thermal transport is that a temperature gradient does 
not correspond to an external ’’mechanical” force like the 
one originating from an electric potential. To bypass this 
problem, time-dependent ’’gravitational potentials” can 
be introducecP^ as source fields in the microscopic ac¬ 
tion. The heat density-heat density correlation function 
can be found by a variation of the action with respect to 
these source fields. Knowledge of the correlation function 
allows to determine the thermal conductivity.^^^ 

An unpleasant difference of the gravitational poten¬ 
tials with respect to, for example, electromagnetic po¬ 
tentials, is that the gravitational potentials couple to 
all terms constituting the Hamiltonian density. This in¬ 
cludes, in particular, the interaction part. Furthermore, 
in the presence of impurities, the gravitational potentials 
also couple to the disorder part of the Hamiltonian. In 
Ref. [51 the latter problem has been overcome by a special 
diagrammatic procedure. Recently, we showed how the 
use of the gravitational potentials can be merged with 
the NLcrM formalism, and performed a renormalization 
group (RG) analysis for the thermal conductivity of a dis¬ 
ordered Fermi liquid system with short-range interaction 
potentials.!^ The RG procedure covers the interval of 
energies with the elastic scattering rate 1 /t as the upper 
cutoff and the temperature T as the lower one (T ^ 1 /t) . 

Combined measurements of thermal and electric con- 
ductivites are often employed in order to assess the ap¬ 
plicability of the quasiparticle description.!^!!^] xhe anal¬ 
ysis of Refs. I4I5I revealed that for the two-dimensional 


such as the density-density or spin density-spin 


disordered system with short range interactions the 
Wiedemann-Franz lavS^ (WFL) holds even in the pres¬ 
ence of quantum corrections caused by the interplay 
of diffusion modes and the electron-electron interaction. 
Generally speaking, the WFL should not be considered 
as a strict law outside the realm of single-particle physics. 
This is already evident from the very fact that the poten¬ 
tial used for calculating the electric conductivity couples 
to the particle density only, while the gravitational po¬ 
tential probes the entire Hamiltonian density. Still, the 
RG analysis shows that at least for the leading logarith¬ 
mic corrections in a two-dimensional system with short 
range interactions, the WFL is obeyed. 

In this paper, we present a perturbative analysis of log¬ 
arithmic corrections to the heat-density heat-density cor¬ 
relation function in a two-dimensional electron gas, i.e., 
in a system with long-range Goulomb interaction. Since 
the effects of the Goulomb interaction in the RG-interval 
of energies have already been studied in perturbation 
theorjEl and are very similar to the case of the short- 
range interaction,]^ we will focus our attention on the 
sub-temperature energy range, which is beyond the scope 
of the RG analysis. The main difference between the RG- 
interval and sub-temperature energy range is that, while 
the transitions described by the standard RG procedure 
are virtual, the sub-thermal range deals with the on-shell 
scattering. For the analysis of the logarithmic corrections 
to electric conductivity, the sub-thermal processes can 
usually be neglected. Thermal conductivity constitutes 
an important exception. Here, the scattering processes 
induced by the long-range Coulomb interaction yield log¬ 
arithmic corrections which, in principle, may compete 
with the RG-corrections. The corrections caused by the 
on-shell scattering, in contrast to those of the RG origin, 
violate the WFL. In this manuscript, we identify the rel- 
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evant diagrams and find the correction to the WFL. We 
show how the terms violating the WFL become compati¬ 
ble with the general form of the heat density-heat density 
correlation function. We thereby demonstrate the con¬ 
sistency of our results with the general scheme for the 
calculation of a correlation function of the density of a 
conserved quantity, which in our case is the energy. 

Our study differs from previous related worlP^^E^lt^ in 
several respects. The heat-density heat-density correla¬ 
tion function was studied before in Ref. [SI However, loga¬ 
rithmic corrections originating from the sub-temperature 
regime were not taken into account in this work. Other 
studies of thermal conductivity available in the literature 
can be divided into Kubo-type linear response calcula¬ 
tions b ased o n the heat-current heat-current co rrelation 
functiorP^E^ and kinetic equation approaches! ^ * ^^ * ^° -*^ 
Our final result for the thermal conductivity of the sys¬ 
tem with Coulomb interaction coincides with the one 
stated in Refs. I8ll81 - f2ni While the mentioned works ar¬ 
rived at the same final result, they did not agree on the 
definition of the heat density and of the associated heat 
current, a question of principle importance for the calcu¬ 
lation of the thermal conductivity. We will devote special 
attention to this point. 

The paper is organized as follows. In Sec. |n]we state 
general properties of the heat-density heat-density corre¬ 
lation function as well as its relation to the quantity of 
our interest, the thermal conductivity. We also introduce 
gravitational potentials as source fields in the action. In 
Sec. [ITT] we define the heat density for the electron system 
with Coulomb interaction and present the NLcrM in the 
presence of the gravitational potentials. This model will 
serve as a starting point for the calculation of the heat- 
density heat-density correlation function. The Coulomb 
problem has a peculiar feature: While we are interested 
in heat transport in a two-dimensional electron system, 
the natural definition of a local conservation law con¬ 
necting heat density and heat current requires a three- 
dimensional setting. The reason is that a part of the en¬ 
ergy of the system is stored in the electromagnetic field, 
and this field is not restricted to the two-dimensional 
plane. In order to define transport of heat in two dimen¬ 
sions, we devise a specific projection procedure. Special 
care has to be taken already on the level of the defi¬ 
nition of the three-dimensional conservation law. The 
principle of gauge invariance plays a pivotal role in un¬ 
ambiguously identifying the heat density and heat cur¬ 
rent. In the present context, this aspect was stressed in 
Appendix B of Ref. [501 We illuminate this point fur¬ 
ther in Appendix]^ where we stress the connection with 
the field theoretic const ructio n of the Belinfante energy- 
momentum tensor .1^21211 In Sec. llVl we collect basic formu¬ 
las required for the calculation of the dynamical part of 
the heat-density heat-density correlation function. Fur¬ 
ther on, in Sec. |V]we discuss the structure of the static 
and dynamic parts of the heat-density heat-density cor¬ 
relation function. Special attention is devoted to the 
consistency with the constraint imposed by the energy 


conservation law. Finally, in Sec. VI we introduce the di¬ 
agrammatic representation and present the general anal¬ 
ysis of logarithmic corrections for the heat density-heat 
density correlation function. In particular, we clarify how 
corrections from the sub-temperature interval of energies 
(caused by on-shell scattering processes) enter the heat 
density correlation function and modify the thermal con¬ 
ductivity. The full list of logarithmic contributions of var¬ 
ious kind is given in Appendixj^ Throughout Secs, rvfvl 


as well as in Appendix we systematically compare the 
heat density-heat density correlation function with the 
well-studied example of the density-density correlation 
function in order to stress differences and similarities. 


II. GENERALITIES: THERMAL 
CONDUCTIVITY AND THE HEAT-DENSITY 
CORRELATION FUNCTION 


In this work, we use the Keldysh techniqu^^lHlZI^ which 
allows us to calculate the correlation function directly in 
real time. The action is defined on the Keldysh time- 
contour C consisting of forward (-I-) and backward (—) 
branches. We start our considerations with the action 

= y dt , (I) 

which contains the heat density k explicitly. The heat 
density is defined as k = h — /in, where h and n are 
the hamiltonian density and particle density and /i is the 
chemical potential. Further, i/; = 
are Grassmann fields with two spin components. 

We wish to calculate the retarded heat density correla¬ 
tion function Xkk{xi,X 2 ) = -iO{ti-t 2 ){[k{xi),k{x 2 )\)Tj 
where x = (r,f), k = ft. — /in is the heat density opera¬ 
tor and the angular brackets denote thermal averaging. 
The definition of the heat density operator reflects the 
fact that we study the propagation of heat under the 
condition when mechanical work (e.g. the radiation of 
acoustic waves) can be neglected. For the calculation 
we define the classical {cl) and quantum components (</) 
of the heat density symmetrized over the two branches 
of the Keldysh contour, = \{k+ ± k-)^ and 

write the retarded correlation function as Xkk{xi,X 2 ) = 
—2i {kci{xi)kq{x 2 )), where the averaging is with respect 
to the action Sk¬ 
in order to generate the heat density correlation func¬ 
tions, we add the source term 


Srj = 2 [r] 2 {x)kci{x) + r]i{x)kq{x)]. 

J X 

to the action. Then, one can find Xkk as 


Xkk{xi,X2) = - 


S^Z 


2 6r]2{xi)Srii{x2) 


1)2 =1)1=0 


( 2 ) 

( 3 ) 


where Z = J '(/)] exp{iSk+iSq) is the partition func¬ 
tion and ip and generalize the definition of ip and ip''' to 
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the Keldysh space. The so-called gravitational potentials 
ryi and 772 in allow us to formulate a linear response 
theory for the heat transport. 

To find the thermal conductivity, it remains to estab¬ 
lish a connection between the response to the gravita¬ 
tional potential and the response to a temperature vari¬ 
ation 5T. As argued by LuttingeilS (see also Refs. H]), 
the responses to 6T and Try may be identified. Namely, 
when considering the response to the gravitational po¬ 
tential, which substitutes the temperature gradient, one 
should replace ry(q, w) —hT(q, a;)/T. For the purpose 
of finding the thermal conductivity from the dynami¬ 
cal heat density-heat density correlation function, it is 
important that the limit q ^ 0 should be taken before 
w —> 0. The heat current jfc(q, w) may be found as a 
response to the gravitation potential ry(q, w) and, fur¬ 
thermore, the static thermal conductivity k will be de¬ 
fined as the real part of the coefhcient relating the heat 
current and —VT using the discussed relation between 
1 ] and ST. Eventually, the thermal conductivity k has 
to be extracted from the disorder-averaged heat density- 
heat density correlation function. As usual, transla¬ 
tional invariance results from the averaging over disorder: 
{Xkkixi,X2))dis = Xkk(.xi-X2). Now, One may introduce 
the Fourier transform of the correlation function, and ob¬ 
tain K as follow^ 


K = -lim lim 

T u ;—^0 yq—^0 


UJ 


ImXfe/c(q,w) 


( 4 ) 


The calculation of the thermal conductivity in this paper 
will be based on this formula. 

The correlation function Xkk obeys the following two 
important relations 


X/cfc(q = 0,a;-)► 0) = 0, (5) 

X/cfc(q 0, w = 0) =-c^T. (6) 

Eq. (§ reflects the conservation laws of energy and par¬ 
ticle number, while Eq. ([^ relates the static part of the 
correlation function to the specific heat per unit vol¬ 
ume at constant chemical potential. 


III. HEAT DENSITY AND COULOMB 
INTERACTION IN THE EXTENDED NLcrM 

The definition of the heat density and the associated 
heat current for the electron gas has been at the center 
of a controversial discussion in recent works on thermal 
transport. Since this question is of fundamental impor¬ 
tance for the calculation of the thermal conductivity, we 
will devote special attention to it. In view of the rela¬ 
tion k = h — /in, and since the expressions for the parti¬ 
cle density and the particle current are well-known, the 
mentioned discussion focuses around the definition of the 
energy (or hamiltonian) density and the energy current. 
At first sight, the answer seems straightforward, as one 
can construct the energy-momentum tensor (EMT) for 


the system of interacting electrons in a canonical way. 
Knowledge of the EMT allows to read off the continuity 
equation relating the energy density and energy current. 
Two problems arise in this context: 

1) The continuity equation obtained from the EMT 
relates a three-dimensional energy density to a three- 
dimensional energy current, while the problem of thermal 
transport for the two-dimensional electron gas requires 
knowledge of two-dimensional densities and currents. 

2) The canonical EMT is not gauge invariant. 

The first point will be addressed in this section, where 
we suggest a simple procedure to project the three- 
dimensional quantities onto the plane. The second 
point, the problem of gauge-invariance, will be addressed 
in Appendix where we remind the reader of the 
field- theore tical construction of the so-called Belinfante 
g]y/[rp|22l23l.^]^jpj^ gauge-invariant expressions 

for the energy density and current. We would like to re¬ 
mark in this context that these expressions have already 
been obtained in Ref. 1201 appendix B without making a 
connection with the Belinfante tensor. 


A. Projection of the electric field onto the charge 
carrying plane 


Before discussing the derivation of the extended NLcrM 
with gravitational potentials, we would like to describe 
the main elements of the projection procedure separately. 
To this end, we will employ the following notation for spa¬ 
tial vectors: r = {x, y)'^ is a 2d vector, r = {x, y, z)'^ is 
a 3d vector, and To = (x, y, 0)^ denotes r embedded into 
the 3d space. We will assume that the 2d electron gas 
(2DEG) is located in the xyy-plane, while the z-direction 
is perpendicular to this plane. We will also use the no¬ 
tation x = (r, t) for a combination of the 2d vector and 
time; for example, the two-dimensional number density 
is n{x) = ipliix- 

The three dimensional hamiltonian density consists of 
a non-interacting and an interacting part h = hg + . 

The transition to the two-dimensional density is straight¬ 
forward for hg. We focus our attention on the interac¬ 
tion p art. I n the Coulomb gauge, it is given as (see for¬ 
mula (AI7| of Appendix 


Amt(ui) = ^ E"(r,t) , 


where E" (r, t) = — V(/3(r, t) and 


T(.L,t) = J dr'' 


;n(r', t) 
r-r'l ■ 


( 7 ) 


( 8 ) 


Here, n{r,t) = n{r,t)S{z), and n denotes the 2d density 
of electrons confined to the 2d plane. Clearly, the field 
Eil is non-zero outside of the 2DEG. In order to obtain a 
two-dimensional energy density, we integrate in the per¬ 
pendicular coordinate z as 


h^ntix) = j 


( 9 ) 
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It is instructive to transform the interaction term 

E” (r, t) = t)VV(r, t) + t). (10) 

Using the Poisson equation —V^(/ 7 (r, t) = 47 ren(r, t)i5(z), 
this decomposition allows us to write 


hint{x) = ^ ^ n(r,t)Uo(i’o - r(,)n(r',t) ( 11 ) 

J^Vo{r-r'^)n{r',t) . 

As a consequence of this integration in z, Eq. 
appears in the second term instead of the original V^. It 
is clear now that interaction part of the Hamiltonian is 
recovered from hint by an integration over the 2 d plane. 


+ 


IOttc^ 


dz 


Hint — 


dv hint{x) 


f n{r,t)Vo{ro-r'^)n{r',t), (12) 

J r.r' 


where Vb(r) = e^/|r| is the familiar Coulomb interaction 
term. 


Returning to Eq. (11), we note that the first term can 


(loosely) be interpreted as a projection of the electric field 
onto the charge it originates from. The second term is 
a correction, for which the point of observation does not 
coincide with the position of the charge. Later, it will 
be shown that the second term on the right hand side of 
Eq. (11) does not contribute to the correlation function 


in the long wavelength limit due to the presence of V^. 
The crucial point here is that the interaction potential 
Vo becomes screened due to the conducting plane. Eor 
the screened potential, unlike for the bare Vq, one can 


neglect the second term in Eq. (11) in the limit of small 
gradients. 


B. Fermionic action with gravitational potentials 


In this section, we prepare the derivation of the NLcrM 
by introducing the gravitational potential into the ac¬ 
tion and further by decoupling the interaction term. Let 
us recall that according to the discussion in the previ¬ 
ous section the full expression for the two-dimensional 
hamiltonian density is h = ho + hint > where hint is given 
in Eqs. ([^ and ([^ and /iq describes propagation of par¬ 
ticles in the presence of disorder 

ha{x) = ■^\/tl:l\/'ipx+Udis{x)n{x). (13) 

Zm 

In order to write the action in the presence of the 
gravitational potentials in a compact form, it is conve¬ 
nient to define a matrix ?)' acting in the space of fields 
tf= (V'+,'|/'-)^ as 


vi + m 0 \ 

0 m-m) 


(14) 


Then, the action S defined on the Keldysh contour can 
be written as 

f)'] = f {idt - [ud^s - m](1 + f]')) ds'if 

J X 

- J + 

JJ drE'^{l+fi')d3^' (15) 

Here, and in the following, we write J) = dt and 
Ix ~ fr f Summation over the spin degrees of freedom 
is implicit. The third Pauli matrix acts in the space 
of forward and backward fields. Erom now on, all matri¬ 
ces acting in the Keldysh space will indicated by a hat. 
Besides, we wrote E' = (E'^,E'i)’^, where 

eE'|(r,t) = -V J dr'Vd(r-r(,)n±(r',t). (16) 


Since our strategy is to project the entire problem onto 
the conducting plane, it will be assumed that rj' = r]'{x) 
does not depend on z. Note that as a result, f] 2 {x) cou¬ 
ples to the two-dimensional heat density. The hamilto¬ 
nian part of this 2d heat density corresponds to the one 
introduced in Eqs. ([^ and (13). , 

Next, the Keldysh rotation can be performed.^^ZHH Xq 
this end, we introduce new fermionic fields 




$ = L&ztf, 



1 -1 

1 1 


■ (17) 


With the help of the two matrices 71 = doj 72 = di in 
Keldysh space, one may form the matrix of gravitational 
potentials fj = 2 hklk- The action after this rota¬ 

tion reads 


J X 

dr E '^(1 + ?7)72E, (18) 

where 

eEfc(r,t) = -V y dr'Vb(r-r(,)4'(r',t)%4'(r',t). (19) 


The last term in Eq. (18) contains four fermionic fields. 
We introduce two real Hubbard-Stratonovich fields 6 * 1 , 2 , 


forming the matrix 9 = J2k=i 2 dfc 7 fc to decouple this 
term. Note that in the case of the Eermi liquid, in 
order to decouple all interaction terms, four Hubbard- 
Stratonovich matrix fields 9^ have to be introduced, 
where the index 1 = 0 — 3 denotes the density and spin 
density interaction channels. Eor the Coulomb problem, 
without account of Fermi liquid-type interactions, only 
the singlet channel, 6 = 0, is involved. For this reason, no 
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index I will be used here. After these transformations, the 
partition function 2 = J exp(iS'[^'^, 0 ]) 

can be written with the use of action 

J X ^ ' 

- [ w^V#(l + 77)V^+ f 0^V-^720, (20) 

J^ ZTTl J 

where 

V~^(x, x') = J dr"dr"'yo"^ (ro - r") 

xV;,(r",r'",t,f')V'(E"'-i-;). (21) 


matrix (r) with respect to its saddle point position 
= At-t', where 

= (^ 0 ‘^-1^ ^ UeasUe, “e = (^ 0 -1 ) 

and = tanh (e/2T) is the fermionic equilibrium distri¬ 
bution function. Here and elsewhere below, 2x2 matri¬ 
ces denoted by the hat symbol act in Keldysh space, with 
the rotation L being already performed. The manifold of 
low-lying gapless excitations is described by rotations 

Q = u o Q o u, Q = U o a^oU, (25) 


Here, fulfils the generalized Poisson equation 


~V((1 -I- 77 (a:))V)H^(r, r', t, t') = 47 re^< 5 (r — r')(5(t — t'). 

( 22 ) 

A useful relation can be obtained for the electron in¬ 
teraction in the action S in Eq. (20 1 


Vn = VoV-^Vo = -{1 + 77 ,^ 0 } + 


1 




VoiW^)Vo. (23) 


For the sake of simplicity, we used a matrix notation 
for the spatial coordinates here. As one can see, the 
above expression reproduces the interaction term given 
in (11). Note that the expression above is not an approx¬ 


imation; there are no higher order terms in 77 . Naturally, 
Vr,=o(a^, a;') = Vo(ro — ~ making contact with 

the theory of the two-dimensional electron liquid in the 
absence of the gravitational potential. Since 772 ( 0 ;) cou¬ 
ples to the heat density, the action given in ( 20 ) taken 


together with the relation (23) reflects the form of the 


hamiltonian density stated by Eqs. (13) and (11). 


C. The extended non-linear sigma model 

In this paper, we concentrate on peculiarities of ther¬ 
mal transport related to the Coulomb interaction. A 
compact description of our approach to the analysis of 
heat transport in a disordered Fermi liquid with short- 
range interactions can be found in Ref. lU while a detailed 
discussion of the NLtrM extended by the gravitational 
potentials was presented in Ref. [51 

As it has been explained in the Introduction, we are 
interested in small energies and long distances. For dis¬ 
tances exceeding the mean free path, the physics is de¬ 
scribed by slow diffusion modes (i.e., modes describing 
density relaxation in the presence of disorder) rather than 
single-particle excitations. Therefore, the fermionic fields 
ijj and tjj^ have to be integrated out. Furthermore, averag¬ 
ing over disorder realizations can be performed assuming 
that disorder is weak in the sense that epT ^ 1 , where ep 
is the Fermi energy and t the transport scattering time. 
Then, the entire physics of the diffusion modes (the so- 
called diffusons) can be encoded in the fluctuations of a 


where U = and {U o U)t,t' = — t'). The 

o-symbol denotes a convolution in time. 

It will be convenient to release the disorder term Udis 
in the action S from the explicit dependence on the grav¬ 
itational potentials. To this purpose, the transformation 
(” A-transformation”) 

A =( 14 - 77 )“^ (26) 


of the fermionic fields was implemented in Refs. I4I5I For 
details of the A-transformation, we refer to these papers. 
As a result of the A-transformation, fj appears through 

the matrix A = 1 — 71771 — 72772 -I- 272771772 -I-- For 

the calculation of the correlation function according to 
Eq. (§, one needs to consider the expansion of A up to 
second order in 7). For the dynamical part of the corre¬ 
lation function, however, only the terms linear in 77 are 
required. 

Starting from the fermionic action displayed in 
Eq. (20), one may apply the A-transformation and subse¬ 
quently follow the traditional route to derive the NLctM 
suitable for description of disordered electrons interact¬ 
ing via the Coulomb interaction. When written in terms 
of deviations of the matrix field Q from its saddle point, 

6Q = Q — A, the model looks as follows 


^="^Tr 


D{WQ)^ + 2i{e,X}SQ 
tT[Xji SQtt ir)]{'j2V^{x, x'))ij 


X tr[A7j- (5Qt/t> (K)] 

+ Tco ( ff{x)^ 2 V{x). 

J X 


(27) 


Here, the tr-symbol includes a trace in Keldysh space, 
an integration over frequencies (when the matrix Q is 
written in frequency space), and a summation over spin 
degrees of freedom; the symbol Tr includes, in addition, 
an integration over coordinates. The first two terms in 
Eq. (27) describe diffusion in the absence of the electron- 


electron interaction; D is the diffusion coefficient; pq is 
the single particle density of states per spin direction. 
The electron-electron interaction acts only in the singlet 
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channel (no Pauli matrices acting in the spin space are 
present) as it should be for the Coulomb interaction. The 
term in the last line describes the contributions to the 
static part of the heat-density heat-density correlation 
function originating from fermionic degrees of freedom 
(i.e., without participation of the diffusion modes); cq = 
271 "^ vqT /3 is the specific heat of electrons. We suppressed 
an additional term that is linear in rj 2 and required only 
for the calculation of the heat density itself. 

The Coulomb interaction entering the action S is stat¬ 
ically screened, V® = (V“^ -I- 2r'oA)“^. This formula is 


symbolical: Both V® and V,, depend on three-dimensional 
spatial coordinates, but screening takes place in the two- 
dimensional plane. Importantly, A appears in the term 
responsible for screening. At zeroth order in 77 , the inter¬ 
action coincides with the statically screened Coulomb 
interaction V®^q = where again 

screening occurs in the plane only. The relation (231 al¬ 


lows us to obtain a regular expansion for in powers of 


V 


= i{l + 7 ), Kf} + (28) 

Since is not singular anymore, one can neglect in 
the second term on the right hand side in the limit of 
small gradients. Thus, owing to screening, the point of 
observation coincides with the position of the charges 
when finding the heat density correlation function of 
a system of conducting electrons confined within a 2d 
plane. (In the case of the bare, i.e., unscreened. Coulomb 
interaction one cannot neglect the second term in Vn 
even in the limit of small gradients.) As a consequence, 
all subsequent considerations involve the effective two- 
dimensional Coulomb interaction with Pq® = 27re^/(|q| -I- 
Ks), where is the inverse screening radius. 


IV. DYNAMICAL PARTS OF THE 
CORRELATION FUNCTIONS - GENERAL 
FORMULAS 


A. Dynamical part of the heat density correlation 
function 


Here, we focus on the dynamical part of the correlation 
function X^|", for which the corresponding diagrams are 
reducible with respect to cutting a single diffuson. The 
starting point for all subsequent calculations will be the 
Keldysh NLcrM action in the presence of the gravitational 
potentials, Eq. ( p7| . Only the Q-dependent part of the 
action ( [27| ) is relevant for the calculation (the last term 
in S can be abandoned). In addition, we may restrict 
ourselves to terms of linear order in t) in the action. This 
allows us, in particular, to use the linear approximation 


for the interaction V® displayed in Eq. (28). To linear 



FIG. 1: The frequency and interaction vertices originating 
from the source terms Sns and S^v as introduced in Eqs. (311 
and (321, respectively. The density vertex arising in connec¬ 
tion with the calculation of the density-density correlation 
function will be drawn in the same way as the frequency ver¬ 
tex. 


order in 7), the Q-dependent part of the action reads 
Sun = ^Tr [d(VQ)2 + 2t{e, 1 - fi}5^ (29) 

- ^ / tr[(l - 7))7,^(r)]7;^Uo®(r - r') 

J rr' 

X tr[7j dQtt (r')]. 

We decompose 


Slin — Srj=0 + ^erj + SjjV ■ ( 30 ) 

with two types of source terms in the action. The first 
one is already present in the noninteracting theory 

=^Tr[{e,77}^]. (31) 

The other source term is specific for the interacting prob¬ 
lem 


SjjV = 




tr[7?7idQtt(r)] 


(32) 


' rr' .t 


X - r')tr[7jA^(r')]. 


As will become clear below, the existence of this vertex is 
of crucial importance for the internal consistency of the 
theory, in particular with respect to the conservation of 
energy. 


The two source terms of Eq. (31) and Eq. (321 give rise 


to two vertices in the diagrammatic representation, which 
we will refer to as the frequency vertex and the interac¬ 
tion vertex, respectively. They are displayed in Fig. 
One can further distinguish between vertices originating 
from a differentiation with respect to 772 and 771 . For the 
sake of definiteness, we will draw the vertices related to 
772 on the left hand side and those related to 771 on the 
right hand side of a diagram. 

Correspondingly, for finding the dynamical part, we 
need to calculate 




3-iii(ei-E2)-|-it2(e3-e4) 


(33) 


X {ei2trYl25Qe-,e2 (i’i)]e34tr[7i(5Qe^£3 (r2)]). 
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where = (e^ + ej)/ 2 , together with the term 


X^^{xi,X2) (34) 

= f e"*‘i(®i"®")(ei2tr[725Q^(ri)] 

° Jr^,€i 

X tr[7i7,5Qt^(r2)]72^V(f (r2 - r3)tr[7j-(5Qt,t, (rg)])^ 


and the analogous term a; 2 ). We introduced the 

notation /^ = / 77 - The index r in these formulas in¬ 
dicates that only those contributions should be selected 
that are reducible with respect to a single diffuson. Av¬ 
eraging (...) is with respect to the action S'^=o- A term 
with two interaction vertices exists, XvVj but is not writ¬ 
ten because it does not contribute to the dynamical part 
of the correlation function in the one-loop approximation, 
but only to the static part. 


B. The dynamical part of the density-density 
correlation function 

It is instructive to compare the calculation of the heat 
density-heat density correlation function to that of the 
density-density correlation function in the same formal¬ 
ism. This correlation function can generated from the 
source term = Tri'QTrlipQ] by differentiation 

Xnn{xi,X 2 ) = , -C , (35) 

2 5ip2{Xi)5Lpi{x2) 

in analogy to Eq. ([^ . From this formula, one obtains the 
expression 

X'}S^{xi,X2) = -|(7ri/o)^/^. e-*‘A£i-£2)-H*t2(e3-e4) 

X (tr[72 dQeie, (ri)]tr[7i(5Q 

£4^3 .(r2)])r. (36) 

In contrast to the heat transport, only a single vertex 
exists, the density vertex. For this vertex, we will use 
the same graphical representation as for the frequency 
vertex, i.e., the one displayed in the left part of Fig. 

For the dynamical part of the density-density corre¬ 
lation function rescattering on the short-range part of 
the electron-electron interaction is allowed, while for the 
heat density correlation function this is impossible. In 
spite of these differences, within the RG-interval of ener¬ 
gies the WFL holds; see Refs. HHSl This indicates that 
there are non-trivial but robust connections between the 
interaction vertices and the interaction amplitudes that 
are fulfilled for the Fermi liquid and remain valid even 
during the course of the RG transformations. 


C. Perturbation theory and dynamical screening 

In the calculation of the correlation functions, an ex¬ 
pansion of Q in deviations from the metallic saddle point 
(J 3 is needed. For the sake of definiteness, we choose 



ip) 


FIG. 2: (a) Diagrammatic representation of the frequency ver¬ 
tex in the sigma model. The first three terms of the expansion 
in P-modes are displayed, (b) Diffuson. 


the exponential parametrization U = exp(—P/2) with 
{( 73 , P} = 0, so that Q = (j 3 exp(P). Fortunately, an ex¬ 
pansion to low orders in the generator P is sufficient for 
our calculation. As an example, the diagrammatic rep¬ 
resentation of the frequency vertex in the sigma model 
is illustrated in Fig. (a); the interaction vertex is rep¬ 
resented analogously. When the expansion of Q in P is 
implemented in Sri=o, this gives among other terms rise 
to the quadratic action Sq of the noninteracting theory 


= -*^Tr 


D{Vpf - 2feV3p2 


(37) 


The action Sq describes the propagation of diffusion 
modes, so-called ’’diffusons” [see Fig. i(b)l 


T> — _ 

o • ' 

Dq^ — iLu 


(38) 


Note that V is the retarded diffuson, the advanced diffu¬ 
son will be denoted as V and is related to the retarded 
diffusion as In the perturbative calcu¬ 


lations involving diffusion modes, Sq given by Eq. (37) 


serves as a starting action. The necessary contractions 
for the Gaussian averaging can be performed with the 
help of the contraction rules that we formulate in Ap¬ 
pendix 

A distinctive feature of the sigma model for interacting 
systems is that upon expansion of Sri=o in deviations from 
the saddle point the interaction potential contributes to 
the quadratic form in the P-modes. This allows one to 
incorporate Fermi liquid effects into the propagation of 
diffusion modes in an automatic way. These effects can 
be interpreted as rescattering of diffusons by the electron- 
electron interaction, or alternatively as a modification of 
the interaction amplitudes by diffusons. An example of 
such a processes is presented in Fig. which represents 
the dynamic part of the polarization operator. [A discus¬ 
sion of the diffuson propagator modified by the electron 
interactions in the Keldysh formalism can be found in 
Sec. Ill B of Ref. [551] In the present context it is more 
convenient to delegate the output of the resummation 
to the interaction itself. As a result, the dynamically 
screened Goulomb interaction should be used instead of 
the statically screened V)f. 



























In the Keldysh formalism, the dynamically screened 
interaction acquires a non-trivial matrix structure in 
Keldysh space: 


= 

k,i/ 


0 


= 




(39) 


In this formula, Vb(k) = 27re^/|k| is the effective 
two-dimensional Coulomb interaction and V^(k.,i') = 
2 i/QD'k^/{Dli^ — iv) is the retarded polarization opera¬ 
tor. The advanced and Keldysh components of V are de¬ 
fined as and - V^^), where 

Bi, = coth{i//2T) is the bosonic distribution function. 


V. STRUCTURE OF THE CORRELATION 
FUNCTIONS 


Before turning to the calculation based on the spe¬ 
cific formalism used in this paper, it is instructive to dis¬ 
cuss the general structure of the heat density-heat den¬ 
sity correlation function Xfcfe(q)W). In particular, we are 
interested in the constraint given by Eq. ([^, which is a 
consequence of the fact that x/cfc(q, uj) describes the prop¬ 
agation of the heat density under the condition when the 
entropy is a conserved quantity. 

We are interested in the singular behavior of Xkki<i,oj) 
which depends on the order of taking the limits q ^ 0 and 
w —> 0. We will assume in this section that all intermedi¬ 
ate integrations have already been performed, and, corre¬ 
spondingly, all corrections arising from the RG-interval 
and sub-temperature energy range have been absorbed 
into the constants which determine the correlation func¬ 
tion. In other words, we will discuss the ” ultimate” stage 
when everything that does not depend singularly on q 
and w can be substituted by a constant. The remain¬ 
ing singular behavior originates from the diffusion prop¬ 
agation of electron-hole pairs, which for free electrons is 
described by the propagator 


I?(q,a;) 


I 

Dq^ — iu! 


(40) 


In the presence of the electron interaction, this propaga¬ 
tor has to be modified as will be described below. 

In order to allow for a direct comparison with the 
density-density correlation function, Xkk will be struc¬ 
tured in the same way as Xnn- [A discussion of the 
density-density correlation function can be found in 
Refs. E5H551 The heat density-heat density correlation 
function has been analyzed in Ref. IHl However, the scat¬ 
tering processes, which are the center of our interest here, 
have not been considered so far for Xkk-] In both cases, 
the correlation function can be split into static and dy¬ 
namical parts. As we have already mentioned in Sec. [TTj 
see Eqs. ^ and Q, the static parts are related to the 
corresponding thermodynamic quantities: the compress¬ 
ibility in the case of Xnn, £^nd specific heat in the case of 
Xkk- The dynamical parts should cancel the static ones in 



FIG. 3: The structure of the dynamical part of the density- 
density correlation function in accordance with Eq. (43l. 


the limit q = 0, a; —>■ 0, which is the way the conservation 
laws for particle number and energy manifest themselves. 
Our goal is to demonstrate how this works for Xkk- We 
start, however, with Xnn for which this procedure is well 
established. 


1. The density-density eorrelation funetion 


The density-density correlation function can be split 
into a static and a dynamical part 

X„„(q, W) = Xnn + Xnr(q>w), (41) 

where the static part is defined as Xnn = Xrm(q —t 0 , w = 
0). Quite generally, the static and dynamical part can be 
further decomposed as follows 


X 


st 

nn 


■ydyn 

Ann 


(q,w) 


-2vqX^ 

- 2 ^ 0 ( 7 .")^ 


^(q, w) -I- iVpUj 


(42) 

(43) 


Let us discuss the parameters appearing in the above 
expressions. As is well known, the static part of the cor¬ 
relation function is related to the compressibility Xnn — 
—dnjdii. Therefore, 


7. 


p — 


1 dn 
2vq dyi 


(44) 


The structure of the dynamical part of the correlation 
function is displayed in Fig. The vertex corrections 
for the two scalar vertices are denoted by 7 ^; fp is the 
short range part of the singlet interaction amplitude. 
This means, in particular, that the long-range part of 
the Coulomb interaction is not included in Xnn- The dif¬ 
fusion propagator modified by the electron interaction, 
the diffuson is defined as 


T’{(q,a;) 


Dq^ — izoj' 


(45) 


It incorporates the frequency renormalization z, intro¬ 
duced in Ref. 1301 and the wave-function renormalization 
Using the relation Xnn(q = 0,w —>■ 0) = 0, a di¬ 
rect consequence of particle number conservation, one 
deduces the following condition 


Zl 




(46) 
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where we use the notation zi for the combination = 
2 ; — p. Adding the static and the dynamical part, one 
finds in view of the Eqs. (46) and Eq. (44) that 


Xrm(q,w) 


dn Dq^ 
dfi Dq^ — iziuj 


(47) 


The wave-function renormalization appears explicitly 
only in the diagramma tic ap proach. In the cr-model ap¬ 
proach to the problem, l2Slini one deals directly with the 
effective amplitude Ep = ^^fp and the effective vertex 
correction 7 ^ = ^ 7 ^. Electric conductivity can be found 
from the relation 


(T = — lim lim 

Cl)— >-0 q—>-0 


f7lniX^„(q, w) 


(48) 


from which one deduces a = 2i'Qe^D. 

Let us now turn to a specific model in which only the 
long-range part of the Coulomb interaction is accounted 
for. [Screening will be included, while the short-range 
Fermi liquid amplitudes both in the singlet and in the 
triplet channel are ignored. By contrast, the short-range 
amplitudes generated by the interplay of the Coulomb 
interaction and disorder have to be included.] It means 
that in the absence of the Fermi-liquid corrections, at 
the initial scale of the RG-integration, i.e., at 1/r, one 
has 7 ^’ = 7 < = 7 ^ = -21 = = 2 = L as well as 

Fp = Fp = 0. From the previous analysis,l2®23l jt is well 
understood that the compressibility Xnn — —dn/dfi does 
not acquire quantum corrections. In the explicit calcula¬ 
tions presented below, we will take the relation 7 ^' = 1 
for the static part for granted, and only analyze the dy¬ 
namical part of the correlation function. The goal will 
be to identify the corrections 6 ^^, 6D, 6z, i^Fp, in the 
expression 


^dyn 

Ann 


(q,w) 


fw(l -I- 6 ^“^ + 2^7^) 

Ziyf) -=— 

{D + SD)q^ — iuj{l + 5z — ^Fp) 


« — 2i'Qiuj'D 

— 2voioj [S^‘^{Dq^ — iuj) — 6Dq^ + iuj{Sz — (5fp)] 

- 2uoiuj{25x^)V. (49) 


In particular, we need to check the relation 5zi = 0, 
which implies 

= 5(e'fp) = 6 fp. (50) 

The first equality follows directly from the definition of 
zi , while the second equality is a consequence of the fact 
that initially Fp = 0 and = 1. A second constraint 
reads 


26 jP = 6^^ + 2 , 57 ^ = 0 . 


(51) 


This constraint follows from Eq. (46) (a consequence of 
particle number conservation) under the condition that 
Szi = 6x^ = 0 . 

Clearly, a full diagrammatic analysis of Xnn for the 
disordered electron liquid requires further steps. Details 
can be found, for example, in Ref. 1311 


2. The heat density-heat density correlation function 

In this section, we will discuss the structure of the heat 
density-heat density correlation function. The discussion 
will be organized in the same way as for the density- 
density correlation function. 

The correlation function can be split into static and 
dynamical parts. The static part is defined as Xkk — 
Xfcfe(q —>■ 0, w = 0), and we write 

Xfcfe(q, w) = xfk + w). (52) 

The static and dynamical parts take the following struc¬ 
ture 



where the propagator 
f>^(q,a; 


-coTxl: 

-coT{xlf[iuj'b^{q,uj)], 

~e 

Dq^ — izuj 


(53) 

(54) 


(55) 


depends on the constants z and and D which have to 
be found during the process of calculation. The relation 
between these quantities and z, and D introduced for 
the density-density correlation function will be clarified 
later. Further, 7 ^ is a correction to the frequency ver¬ 
tex. The static part of the correlation function, as follows 
from Eq. (§, is determined by the specific heat c of the 
electronic system, which is known to acquire quantum 
corrections within the renormalization group interval of 
energies. In Eq. (53), these corrections are absorbed into 
the quantity 



The structure of the dynamical part of the heat 
densi^-heat density correlation function is displayed in 
Fig. W Its main difference from Xnn" is that for xS" 
ladder diagrams with the interaction amplitudes Fp as 
shown in Fig. are not relevant. Therefore, the singu¬ 
larity of this correlation function is determined by the 
denominator of the diffuson propagator without in¬ 
sertions describing rescattering; compare Fig. Elto Fig .0 
The reason underlying this observation can be under¬ 
stood easily. An insertion of the static amplitude Fp de¬ 
couples the frequency integrations on the left and right 
hand side of the diagram. As a consequence, a frequency 
integral of the type £{J'e+uj/ 2 —J^s-uj/ 2 ) = 0 arises from 
the vertex related to the classical component of the grav¬ 
itational potential, and diagrams with Fp-insertion do 
not contribute to xtT^- important to note, however, 
that this simple observation does not imply the absence 
of vertical diagrams in general, as will be discussed in 
detail in the next section. 

Energy conservation imposes a constraint on the cor¬ 
relation function Xkk, which is encoded in the relation 
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FIG. 4: The structure of the dynamical part of the heat 
density-heat density correlation function in accordance with 
Eq. (pi. 


Xfcfc(q = 0,a; 0) = 0, compare Eq. This, in turn, 

imposes the following constraint on the parameters en¬ 
tering the static and dynamical correlation functions 

= (57) 


As has been discussed in Refs. [Ml the quantity z which 
describes renormalization of the frequency term in action 
S', and in this way enters the propagator of diffusons, is 
directly related to the specihc heat, c = zcq, so that 
= z. Using this information as an input, Eq. (57) can 
be also written as zz = ( 7 ^)^, where we defined 7 ^ = ^ 7 ^. 
Within the renormalization group interval of energies, 
this relation degenerates to z = z = 7 ^. Adding the 
static and the dynamical parts, one then finds 



FIG. 5: Diagram for the non-interacting part of the dynamical 

1 , • r 1 - dyn 1 dyn 

correlation functions Xee^o Xnn,o* 


VI. DYNAMICAL CORRELATION FUNCTIONS 
- DIAGRAMMATIC ANALYSIS AND 
LOGARITHMIC CORRECTIONS 


In this Section, we present an analysis of the dynamical 
part of the heat density-heat density correlation function 
in the diffusive limit. The analysis will be based 
on the NLcrM action derived in Sec. Ill Eq. (27). To 
highlight similarities and differences, we contrast the cal¬ 
culation of Xfefc” with that of Xnn" within the same for¬ 
malism. In order to prepare the discussion of the inter¬ 
action corrections, we first summarize the results for the 
non-interacting case. 


A. The non-interacting part of the correlation 
function 


Xfcfe(q.w) 


-llcoT 


Uq 2 


Da^ — izoj 


(58) 


In the absence of interactions, only the frequency- 
frequency correlation function Xee” contributes to Xfcfe^j 


At the scale I/r, the initial values for the various pa¬ 
rameters of the theory are 7 ^ = 7 ^ = .S = ? = 1 , and 
the propagator of the diffuson is equal to I?(q, w), com¬ 
pare Eqs. (45) and (55) with Eq. (40). Coming back to 


the dynamical part, Eq. ( |54| ), we therefore expect that 
a perturbative calculation of the dynamical part of the 
correlation function will result in an expression of the 
following form 


xS"(q>w) ~ -cqT 


iuj{l + + 25'y^) 

{D + 6 D)q^ — z(I -I- Sz)uj 

-coTiioidi'^ + 2dr^)V^,c. 

-CoTiui[-6Dq^ + icoSzpl^^. (59) 

To check consistency of the sum of dynamical and static 
parts of Xkki^^^) with the conservation laws, one should 
make certain, in view of Eq. (57), that the relation 6 ^'^ + 
2 (57^ — 6z = Sz indeed holds. 

In the next section, the analysis of the logarithmic cor¬ 
rections to Xkk as well as Xnn is presented. In particular, 
in the following sections. Sec. jVI Aj and Sec. |VIB| the 
structure of the different terms is discussed together with 
their diagrammatic representation, while in Sec. jVI Cj 
logarithmic corrections arising from the RG and sub¬ 
temperature intervals are described in detail. In Ap¬ 
pendix 0 a comprehensive list of different contribution 
is given. 


xS.o(a^i- 2 ; 2 ) = 


^-iti(ei-e2)+it2{e3-e4) 


X (ei2tr[72(T3T^(ri)]e34tr[7i(T3^^(r2)])o. (60) 

The corresponding diagram is displayed in Fig. With 
the help of the contraction rules (B2) or (|B3|) on finds 






(61) 


Here, we introduced the window function 


^e,u; — s+uj{2 ^ e— a;/2- 


(62) 


The appearance of the window function is characteristic 
for the dynamical part of the correlation function. For 
T —>■ 0, it allows frequencies e to lie in the interval (e — 
a;/2, e-|-a;/2); at finite temperature this range broadens. 
Still, upon integration in £, the function gives rise 
to the factor of w. Returning to the calculation of Xee oj 
after expansion in w and with the help of the relation 
= ttT'^ I one obtains 

XeS(q,w) = -coTiujVq^^, (63) 


where we remind that cq = 2tt'^i'qT/3 is the specific heat 
in the absence of quantum corrections. 

In complete analogy, one can calculate the dynami¬ 
cal part of the density-density correlation function in the 
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non-interacting limit, 

= -^iTTVo'Dq^^ (64) 

Using the relation tt = to, one finds 

(65) 

The diagrammatic representation for Xnn"o coincides 

with the one for xts'o^ compare Fig.j^ This is the origin 
of the WFL in the case of non-interacting electrons. 


B. Interaction corrections: Diagrams 

We now turn to the explicit calculation of quantum 
corrections to the correlation functions originating from 
the combined effect of the long-ranged Coulomb interac¬ 
tion and disorder. The calculation is performed using an 
expansion in deviations SQ from the saddle point, and 
applying subsequently the contraction rules formulated 
in Appendix]^ Diagrams are presented only for illustra¬ 
tion. A detailed account of the calculation is presented in 
Appendix Here, we will highlight the most important 
diagrams and summarize the results. For comparison, we 
present the information for Xkk in parallel with Xnn- 
We will group the relevant diagrams for the calculation 
of the correlations functions into five classes. When we 
draw the diagrams, we leave out additional partner dia¬ 
grams that can be obtained by simple symmetrization of 
those already displayed. 

1. Horizontal diagrams: These diagrams contain a 
horizontal interaction line and give rise to correc¬ 
tions to the diffusion propagator. Vertex correc¬ 
tions with horizontal interaction lines will be con¬ 
sidered separately. The horizontal diagrams are 
displayed in Fig. The corresponding corrections 
will be labeled as xtTi Xnnp- 

2. Vertical diagram: The diagram with vertical inter¬ 
action line relevant for our calculation is displayed 
in Fig. [Tj It results in corrections to the diffu¬ 
sion propagator. Vertex corrections with vertical 
interaction lines will be considered separately. The 
vertical diagram leads to the corrections x‘kk ^2 

dyn 
Xnn,2 ■ 

3. Drag diagrams: The drag diagrams contain two 

screened interaction lines and give rise to correc¬ 
tions to the diffusion propagator, see The 

resulting corrections will be labeled as XkT^s 

Vertex corrections of drag-type will be con¬ 
sidered separately. 

4. Regular vertex corrections: In this class, we sum¬ 
marize those vertex corrections that originate from 
the frequency vertex Sr/e- Horizontal and vertical 



FIG. 6: The four horizontal diagrams contributing to xtt\ 
and Xnn”i- Each diagram has a symmetry-related partner 
that is not displayed here bnt accounted for in the analytical 
expressions discussed in the text. 



FIG. 7: The vertical diagram which contributes - together 
with its symmetry related partner - to Xkk "2 XX'X, 2 - 


(regular) vertex corrections are displayed in Fig. 
(regular) vertex corrections of the drag type in 
Fig. The regular vertex corrections will be re- 

ferred to as XkkA XnnA- 


5. Anomalous vertex corrections: The anomalous ver¬ 
tex corrections result from the interaction vertex 
generated by Sr^v- Obviously, they only arise in 
the calculation of the heat density-heat density cor¬ 
relation function. Fig. im shows the diagrams for 
anomalous vertex corrections with a single interac¬ 
tion line. They will be labeled as Xu^ 5 - Fig- 


12 


shows the diagrams for vertex corrections of tne 
drag type, i.e., with two interaction lines. These 
corrections will be labeled as xtHTe- 


C. Analysis of logarithmic corrections 

In this Section we will compare logarithmic corrections 
to the diffusion coefficient and the frequency renormal¬ 
ization arising in the heat density-heat density correla¬ 
tion function with those in the density-density correlation 



FIG. 8: The four drag diagrams which contribute to xtT^a 
and together with their symmetry-related partners. 
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FIG. 9: The regular vertex corrections, xS '4 Xnn^- 
Three more diagrams are obtained by symmetrization. 



FIG. 10: Drag-type diagrams for the regular vertex correc¬ 
tions. An explicit calculation shows that their contribution 
to both and vanishes. 


function. The logarithmic corrections can be classified 
according to the most important frequency and momen¬ 
tum regions in the integrals: 

For (i)-a terms the frequency v transferred by the elec¬ 
tron interaction is larger than the electron frequency |e| 
as well as temperature. The frequency integrations are 
controlled by the combination Te+u ~ ^e-v The trans¬ 
ferred frequency and momentum cover the whole RG- 
interval. Electron-hole pairs excited via the interaction 
are virtual and contributions from the sub-thermal region 
are insignificant. 

For (i)-b terms the frequency transfer is limited by the 
combination In this case, the frequency 

transfer is insignificant but the momentum integration 
covers the whole RG-interval. These terms describe the 
modification of the interaction amplitudes by disorder. 

Finally, there appear new contributions, (ii) terms, 
which are determined by the combination J-e-v 

In this case, the transferred frequency is limited either 
by temperature or by jej. Furthermore, the electron 
interaction enters the integrals via its imaginary part, 
This, together with the fact that the transferred 
frequency is limited either by e or temperature, indicates 
that inelastic processes intervene. The momentum inte¬ 
gration is determined by small momenta. 

Logarithmic integrals appearing in (i)-a and (i)-b terms 
will be denoted as li-integrals; see Sec. VIC2| and Ap- 



FIG. 11: Two anomalous vertex corrections contributing to 
Xfefe" 5 - No analog exists for Xrm ■ Two more diagrams are 
obtained by symmetrization. 




FIG. 12: Drag-type contributions to the anomalous vertex 
corrections xtT^e- No analog exists for Xnn™- Two more dia¬ 
grams are obtained after symmetrization. 


pendix They are well known from the previous RG 
studies of the the disordered electron liquid. In contrast 
to electric transport, the contributions (ii) are specific 
for thermal transport; they are important in the case of 
the Coulomb interaction when is singular. For 

a given frequency iz, most important momenta fulfill the 
inequality \iy\/{DKs) < k < y/\v\jD. In this interval, one 
can approximate the dynamically screened interaction as 


ImK^ 


k.i- 


1 V 


( 66 ) 


Eventually, the bare l/Hk^ singularity gives rise to log¬ 
arithmic corrections. T hese log arithmic integrals will be 
denoted as see Sec. 


VIC 3 


and Appendix The 


dex h emphasizes their importance for heat transport. 

Note that the interval \iy\/{DKs) < k < sJlvl/D is 
also responsible for the double-logarithmic dependence 
of the tunneling density of states as well as other spu¬ 
rious corrections that appear in intermediate stages of 
the RG procedure (compare the integral Ii introduced 
below). In the case of the (ii)-type integrals, however, 
only a single logarithm arises, because allowed frequen¬ 
cies v are small, of the order of the temperature, while 
a double-logarithmic dependence is obtained for an (i)- 
a type integral /i, where the frequency can take large 
values. 

In summary, we encounter two different types of con¬ 
tributions. For the first type, which includes (i)-a and 
(i)-b terms, at least one of the two energies |zz| and 
lies in the RG interval (T, 1/r) giving rise to logarithmic 
integrals Ii. These corrections are well studied for the 
case of the density-density response function, both on a 
diagrammatic level and on the level of the field-theoretic 
NLcrM. Concerning the heat density-heat density correla¬ 
tion function, a diagrammatic study has been presented 
in Ref. [HI while the NLtrM of Ref. I4I5I focused on the RG 
in the disordered Fermi liquid, i.e., in a disordered sys¬ 
tem with short-range Fermi liquid-type corrections. A 
common result of these studies was that the logarithmic 
corrections originating from the RG interval for the heat 
density-heat density correlation function lead to the se¬ 
quence of equalities z = = z = 

The second type of logarithmic corrections, the (ii)- 
terms which originate from sub-thermal frequencies v < 
T, are at the center of our interest here. For these cor¬ 
rections, the imaginary part of the dynamically screened 
interaction is relevant. 

Besides these two types, there are terms that could, 
in principle, introduce a mass into the diffuson. They 
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will be denoted as J-terms, for details see appendix [C| 
Unlike U and which enter the calculation of the dy¬ 
namical part of the heat-density heat-density correlation 
function together with the factors Dq^ or uj only, these 
terms are finite (i.e., they do not vanish) in the limit 
(q, w) —?► 0. The J-terms arise as fragments of individual 
diagrams, but they have to cancel in the overall result 
for the correlation functions of conserved quantities. If 
they would persist, this would lead to a violation of the 
conservation laws. The cancellation of these terms is in¬ 
timately related to the balance between in and out-terms 
in the collision integral integrated over frequencies (and 
for the case of the heat density correlation function also 
weighted with frequency). 


1. The finite J-terms 

We start with the J-terms. For the corrections arising 
in the density-density correlation function one gets: 

Xun,ij = -2j^o*wJi(q,a;)X>2 

=-2i"o*wJ2(q,a;)X>2 (67) 

These formulas are obtained from the expressions given 
in appendix after performing the integration in the 
electronic frequencies e. Upon expansion in Dc^ and 
w, Ji(q, w) = J° -|- one notices that Ji 

and J 2 both contain non-vanishing constant parts and 
singular expansion coefficients Jf’ and J“. It turns out, 
however, that there is a full cancellation between hori¬ 
zontal and vertical diagrams: J 2 (q, w) = — Ji(q, w). In 
particular, the cancellation between and J® = — J® 
ensures that the density-density correlation function re¬ 
mains gapless. It is instructive to interpret the cancella¬ 
tion in the limit (q, w) —?► 0 in the language of kinetics (for 
a more detailed discussion see appendix[D|). It can be seen 
that it is a direct result of condition f^SIcoiii£,x) = 0 
for the linearized collision integral, which ensures the con¬ 
servation of the particle number in a kinetic formulation 
of the problem. 

Next, let us look at the corrections arising in the heat 
density-heat density correlation function: 

= -CoTiujJi{q,u})Vl,^, 
xtl^2j = -CoTiuj{J2{q,U}) + J2)T>q,<^, 
xtZj = -coTtujJ.Vl^. ( 68 ) 

Here, the new terms J 2 and J 3 , are finite and do not 
contain q and w dependent parts. The cancellation of 
the J-terms, which results from the identities J 2 = — J 3 
in addition to J 2 = — Ji, is now somewhat more com¬ 
plicated and involves the drag diagrams. In the present 
case, the cancellation is not guided by the number con¬ 
servation or, equivalently, the absence of a mass of the 
diffuson, but by energy conservation. In the language of 
kinetics it can be seen that the cancellation is a direct 


consequence of the condition ^ e6Icoii{£, a;) = 0 for the 
linearized collision integral, which ensures the conserva¬ 
tion of energy. 

Due to the importance of the identity Ji = — J 2 , and 
a similar relation J 2 = — J 3 for the heat density-heat 
density correlation function, we devote Appendix]^ to a 
more detailed discussion of this point. This discussion 
elucidates the relation between the horizontal, vertical 
and drag-type diagrams in the low-energy interval. 


2. Logarithmic corrections from the RG interval 


Here, we list the logarithmic corrections to Xnn"o = 
—2i'iuj'D originating from the RG interval. They are en¬ 
coded in the logarithmic integrals denoted as U. Detailed 
derivations as well as the definitions of the appearing in¬ 
tegrals li can be found in Appendix [C] 

Xnn”i = -2r'o*w [-2{Dq^ - iu})Ii Dq^Io - iuth] 2?^, 

xtT,2 = - 2 r'o*w {iujl^) 

dyn ^ „ 

Ann,3 

xtT,4 = -2r'o*w/iX>. (69) 

In this list, we suppressed the arguments of V for the 
sake of brevity. Concerning the vertex correction 
the list presents the correction for one individual vertex, 
i.e., the sum of the vertex corrections from the left and 
the right vertex is twice as large. 

For the heat density-heat density correlation function, 
we find the following corrections to xS"o “ — cqT 11012 


XfeM = -CoTioo [-2{Dq^ - iuj)h + Dq^Io - iool,] 
xtk^2 = -coTiw [iw{h - h)] 

dyn _ g 

xS?4 = -coTiujhV, 
xtT ^5 = -coTiu:{-h)V, 

xtT, = 0. (70) 


We would like to stress that those integrals in Eqs. (69) 
and (70), which are denoted by the same names, are not 
only equal but determined by the same expressions. Sim¬ 
ilar to the case of Xnn"j the list cites vertex correction for 
one individual vertex only. 

Due to the presence of the window function « 

one may set the electron frequency |e| « T in the 
expressions determining the integrals A. Then one finds 


h = ^Plog ^ log (71) 

lD=plog^, (72) 

Iz=l2=h=\lD, (73) 
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with p = (A'K'^V[jD)~^. Whereas the main contribution 
for the momentum integral in Ii comes from the inter¬ 
val \v\/{Dks) < k < ^/\v\JD, for the rest of the terms 
relevant momenta are such that > \'>^\- Relevant fre¬ 
quencies are large > T for /i and Id, but < T for 


Iz, h, and h. 

Density-density correlation function: Following the 
logic of Sec. |Vj we arrange the obtained corrections to 
Xnn into the general form consistent with that of a cor¬ 
relation function of a conserved quantity. Then, from 
comparison with the results listed in Eq. (69), and the 
expression Eq. (49) we can find the corrections to the 
various constants characterizing this correlation function. 
First of all, we observe that the wave function renormal¬ 
izations given by the /i-term in and the vertex 

corrections Sx<i given by Xdi^ 4 , cancel out: 


= -2/i (74) 

6Y^=h. (75) 


This ensures the absence of doubly-logarithmic correc¬ 
tions in Xnn- Furthermore, the effect of the frequency 
renormalization, 6 z = —Iz, given by the last term in 
Xnn"i i® cancelled by that of the renormalized screened 
Coulomb interaction, JFp = —Iz, given by Xnn,2- The 
only effective correction which remains after the cancela¬ 
tions is a correction to the diffusion coefficient in Xnn"i ■ 

SD = -DId. (76) 


Thus, we reproduced the following (known) results: 

1. The density-density correlation function in the pres¬ 
ence of quantum corrections (albeit ignoring Fermi-liquid 
type corrections) reads 

D 

X„„(q,w) =- 21/0 7^— 2 ■ > (77) 

where the diffusion of charges is governed by the charge 
diffusion constant = D -\- SD. 

2. Electric conductivity can be found from the relation 


(T = —lim lim 

cj —>-0 q —>-0 


— ImXrin(q,w) 


= 2vQe^Dn. (78) 


As a consequence of Eq. (78), the correction to con¬ 
ductivity is 


5<J 6D 
a D 


-Id = -plog 


1 

TV' 


(79) 


In this way, one recovers the well-known Altshuler- 
Aronov correction to conductivity from the formalism. 
This correction originates from the RG interval of ener¬ 
gies. 

Heat density-heat density correlation function: A com¬ 
parison of the results listed in Eq. (70) and the general 
corrections stated in Eq. (59) leads us to the following 


relations for the corrections originating from the RG in¬ 
terval 


^ 

SD = SD = -DId, 

Sz = Sz = —Iz, 

sr< = h-h- (80) 


As it has been discussed in Sec. |Vj for the consistency 
of Xkk (q, w) with the energy conservation law, the condi¬ 
tion Sf,'^ -I- 2^7^ — Sz = Sz is necessary. This condition is 
fulfilled provided that Iz = h- While Iz and T 5 are a pri¬ 
ori different integrals, they do coincide with logarithmic 
accuracy and the relation stated in Eq. (57) holds. We 
would like to stress that T 5 originates from the anomalous 
vertex correction, which only exists for Xfcfe"^ (and is ab¬ 


sent for xi’^T) and, therefore, the presence of the source 
term SrjV as already mentioned is a very important in¬ 
gredient of the theory. 


Thermal conductivity can be found from the formula 


1 


w 




IcoD. (81) 

z 


In the last equality we used the form of the correlation 
function stated in Eq. (58) as well as the relation c = cqZ 
introduced before. This implies the relation 


K ZCqD 

aT z2vQe^DnT 


zb 

zDn 




(82) 


where Cq = 7 r^/ 3 e^ is the so-called Lorentz number. 

We thus arrive at the following conclusions concerning 
the Wiedemann-Franz law: 


1. If there were no additional corrections from the sub¬ 
temperature interval, then the set of equations listed in 
( [80| ) would immediately lead us to the conclusion that 
the WFL is fulhlled. Indeed, as one can see from the 
second and third relation in (80), the equalities z = z 


and SD = SD hold. Then the WFL remains true even in 
the presence of the quantum corrections originating from 
the RG interval. 


2. In order to obtain a violation of the WFL, the inequal¬ 
ity zDn zD is required to hold. 


3. Logarithmic corrections from the sub-temperature 
interval 


For each diagram, only the corrections from the sub¬ 
temperature energy interval are listed below: 
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dyn _ 

A/cfc,l ~ 

dyn rri ■ 

Xkk ,2 = -cqTiuj 


-{Dq^ - iuj)i^ 
xtZ = -CoTzu; [l^Dq^ - 


,2 


xtTi = -coTiuj-I^V, 

dyn __ „ 

A/cfc,5 ~ ^ 

xtZ = -coTzuz{-I^)V. 


(83) 


Only corrections from the sub-temperature regime are 
discussed here, therefore we can set 5z = Q when check¬ 
ing the consistency with the energy conservation law for 
Xfcfe(q, w), which reduces to 25^^ — Sz = 0. 

Finally, according to Eq. (81), the correction to ther¬ 
mal conductivity reads 


, r, 1 T ^ DkI 

OK = - log - 1 - log - . 

6 ^Tt 12 ^ T 


4- Violation of the Wiedemann-Franz law 


( 86 ) 


Note with respect to XkCi xtT^e^ that the list cites 
vertex corrections for one individual vertex only. Con¬ 
cerning the integrals , we notice that again we can set 
|e| ~ T and then 

nf-2 

I^=p\og^ (84) 


and I 2 = I 2 = I 4 = 2.1^ = 2/g = . In these integrals, 

relevant momenta are in the interval \v\/{Dks) < k < 

Sec. 


/D and relevant frequencies are small \v\ < T; see 


VIC for a general description of the /f-terms, and 


Appendix [C| for their detailed analysis. 

Since in this manuscript, we study only first order log¬ 
arithmic corrections to the correlation function, the accu¬ 
racy of the calculation is not sufficient to make a definite 
statement about the structure of the correlation function 
as a whole. In particular, unlike for the RG corrections, 
the classification of the corrections in terms of dH, 5z 
and d 7 <| is not unambiguous. This remains so even if 
we assume for the wave function renormalization that it 
is unchanged, 5^ = S^, and take into consideration that 
the specific heat c/cq = z is not affected by the sub¬ 
temperature corrections. There still remains a degree of 
freedom for 7 ^ and z within Eq. (57). The final result, 
of course, will not depend on the choice of presentation 
of the correlation function Xkk- 

Here, we fix the ambiguity following the origin of cor¬ 
rections in Eq. (83). Then, the vertex corrections given 

by Xkk^i and xtlTe cancel in total, 5% = 0. Next, the fre¬ 
quency corrections to 5z originating from the frequency 
terms in and xSis also cancel, so that 5z = 0. 

Thus, this procedure leads us to the following set of sub¬ 
temperature corrections: 


= (5%^ = 0 , 

5z = 0, 

gjjh = 1 /^. (85) 

In the procedure chosen here for fixing parameters, the 
structure of the correlation function (i.e., the vertex cor¬ 
rections and frequency renormalization) are controlled 
by the RG-interval, while the heat diffusion constant be¬ 
sides the corrections from the RG interval acquires a spe¬ 
cial contribution from the sub-temperature energy range, 
SD^^. 


From the results collected in this section we can draw 
the following conclusions: 

1. The full heat density-heat density correlation func¬ 
tion can be written as 

Xkk{fl,i^) =, (87) 

- luj 


where Dk = {Dn+5D^)/z is the heat diffusion con¬ 
stant, and c = zcq. The form presented in Eq. (87) 
is canonical for a correlation function of a density 
of a consered quantity in the presence of disorder. 


2. Comparing corrections to the heat and electric con¬ 
ductivities 


5k 5a 1 u 

— = — + 

K a 2 

one gets that the Lorenz ratio is enhanced 
CoaT 2 


( 88 ) 


(89) 


with = p\og{DKl/T) > 0 . 

The positive sign of the correction indicates that for the 
disordered electron gas with long range Coulomb inter¬ 
action heat transport is more effective than the WFL 
suggests. 


VII. CONCLUSION 

We conducted an analysis of the heat-density heat- 
density correlation function in order to obtain quan¬ 
tum corrections to the thermal conductivity of the dis¬ 
ordered electron gas. Our analysis focused on the role 
of the long-range Coulomb interaction in the diffusive 
limit and combined effects originating from different en¬ 
ergy scales. RG-type corrections arise from the energy 
interval from the elastic scattering rate down to temper¬ 
ature (I/r ^ T) and do not violate the WFL. Additional 
corrections to the thermal conductivity originate from 
the sub-temperature energy range. These corrections do 
not have an analog for electric conductivity and there¬ 
fore violate the WFL. The resulting thermal conductiv¬ 
ity exceeds the prediction of the WFL. In this sense, heat 
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transport is more effective than charge transport. This 
result should be contrasted with the case of the disor¬ 
dered Fermi liquid, a model system with short range in¬ 
teractions. In the disordered Fermi liquid the WFL law 
is obeyed. 

As for the origin of the WFL violating corrections we 
would like to stress that (i) the scattering processes in 
question involve on-shell energies ~ T as well as the 
sub-temperature energy range and (ii) the relevant log¬ 
arithmic integrals contain the imaginary part of the dy¬ 
namically screened Coulomb interaction Imld^(k, i/). We 
therefore conclude that inelastic scattering is at the ori¬ 
gin of the violation of the WFL. We checked that the 
correction to the heat diffusion coefficient caused by the 
long range Coulomb interaction is not modified by Fermi 
liquid interaction amplitudes. We thereby expect that 
the answer obtained for the correction to thermal con¬ 
ductivity Sk presented in Eq. ( 86 ) is final. 

We studied thermal conductivity in a situation where 
mechanical work (e.g., radiation of acoustic waves) can 
be neglected. If one additionally takes the conservation 
of particle number into consideration, this implies that 
heat transport is to a large extent governed by energy 
conservation. Special care has been taken regarding the 
definition of the energy density in the presence of the 
long-range Coulomb interaction. As the energy density 
depends on the electric field, the natural definition of a lo¬ 
cal energy conservation law requires a three-dimensional 
setting. For finding the three-dimensional energy density 
we used the field-theoretic construction of the energy- 
momentum tensor in combination with the principle of 
gauge invariance, which was used to lift the remaining 
ambiguity. These considerations n atura lly led to the Be- 
linfante energy-momentum tensor pinally, in order 
to define an effective two-dimensional energy density, we 
employed a projection onto the plane. 

It is instructive to contrast the thermal conductiv¬ 
ity in the diffusive limit studied in this paper with 
known results in the clean electron gas with Coulomb 
interaction^^ or in the ballistic limit.l^ In the latter 
cases, inelastic scattering processes are responsible for 
a decrease of the thermal conductivity. In contrast, the 
corrections in the diffusive limit lead to an increase of 
the thermal conductivity. The positive sign of the cor¬ 
rection indicates that the incoming scattering processes 
are dominant. Loosely speaking, in the diffusive case 
with long-range Coulomb interaction, electrons can use 
the energy ^ T from a remote region to facilitate heat 
transfer. 
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Appendix A: On the gauge invariant definition of 
the heat density 

The purpose of this appendix is to derive gauge- 
invariant expressions for the heat density and the heat 
current in the presence of the long-range Coulomb inter¬ 
action. To this end, we follow the general proc edure for 
the construction of the Belinfante tensor ,l22l2ll which is 
used for the energy-momentum tensor in electrodynam¬ 
ics. We start with the Lagrangian density (Schrodinger 
field coupled to electromagnetic field) C = Cs + Cemi 
where 

‘-s = \ (Al) 

-^(iV — qA)^* — qA)ijj — qcjitl;* ij} 

2 m 

is the Lagrangian of the non-relativistic Schrodinger field 
Ip with charge q and mass m coupled to the electromag¬ 
netic field A), and 

Cem = (A2) 

Idtt 

is the Lagrangian of the free electromagnetic field.l^The 
potentials (p and A are related to the electric and mag¬ 
netic fields 


E = -V0 - A, B = V X A. (A3) 


Also, = d^A'' — d^A^ is the field strength tensor. 
The relativistic notation is used for convenience only. 
The equations of motion obtained by a variation of the 
action S = f dxC with respect to ip*, ip, A° = cp and A 
give the Schrodinger equation 

idtip = - qA)'^ip + q(pip, (A4) 

and its conjugate, and the Maxwell equations VE = 47rp 
and V X B — 5(E = 47rj, respectively. Here, we defined 
the charge density p = qip*ip and the current density 

j = ^ - qA)ip + ((iV - qA)ip*)ip]. (A5) 


The remaining two Maxwell equations, VB = 0 a nd V x 
E -I- 9tB = 0, are fulfilled automatically through (A3). 


The canonical energy momentum tensor 0^“^ is ob¬ 
tained using the invariance of the action with respect 
to the translation = x'^ + 


= 


dC 


d{d^,ip) 


d''ip- 


dC 


d{df,ip*) 


d^ip* 


dC 


did^A-) 




(A 6 ) 
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where = diag(l, —1, —1, —1). We know that 

= 0. It means that for each i/ we get a local con¬ 
servation law (continuity equation). The conservation 
law related to the energy density is given by 9^0^° = 0 . 
We, therefore, should calculate the energy density 
and the i-th component of the energy current 0®°. One 
finds 




-E9tA-t- 


(A7) 


where we defined 

- qA)'ilj*{-iV - qA)!/;. (A 8 ) 

Zm 

With the help of the Maxwell equations one can rewrite 
this result alternatively as 


0™ = u^ + ^{E^ + B2) + ^V(<))E). (A9) 


Next, we turn to the components of 0*°, for which we 
find 




where 




i 

2 m . 


— qA)ifj — (zV — qA)il;*dt'ip ■ 

(AlO) 


Note that is a gauge-invariant quantity, as a local 
phase change of tp, ip* can be absorbed by A. We 
conclude that the energy density is gauge invari¬ 
ant. As to the current, it can be easily checked that 
the combination jg’* — (/>j* is also gauge invariant, i.e., the 
transformation ip —>■ exp(zx)^ can be compensated by 
A —>■ A -f q~^Vx <p ^ (p — q~^dtX, and so is T*°. 

Note that in the absence of external fields, and neglect¬ 
ing fluctuating magnetic fields, which is a relativistic ef¬ 
fect, we find agreement between the components of 
the canonical energy-momentum tensor and the interme¬ 
diate expressions jf and u? considered in App endix B of 
Ref. [201 as well as between our expressions (A14| and 
(A15) for the gauge-invariant energy density and current 
and the final expressions obtained in Ref. 1201 

Next, we specialize on the Coulomb gauge, the gauge 
used in the main text. It is convenient to decompose 
A = All -f A-'-, with VA^ = 0 and V x AH = 0. The 
Coulomb gauge VA = 0 eliminates the longitudinal de¬ 
grees of freedom AH = 0, so that A = A-*^. The elec¬ 
tric field E = —VAq — dtA, in contrast, has both a 
longitudinal and a transversal part, EH = —VAg and 
E-'- = —dtA-^. Then, Aq is determined by the Poisson 
equation and 


eII 



|r - r'l' 


(AI 6 ) 


Again, with the help of the Maxwell equations, one can 
find the alternative representation 

0" = + (All) 

+ 1-[V X ifB) - dtiEp)]\ 

As is well known, there is a problem with the canoni¬ 
cal energy-momentum tensor; it is neither symmetric nor 
gauge invariant. This is already obvious from the terms 
E9(A and (B x dtA) in the expressions for 0°° and 0°*, 
respectively. Since the densities and currents are not 
defined uniquely, one can add a four-divergence to the 
energy-momentum tensor as 

Tfiu = -f l9‘^Xcr/j;y, (-^12) 

where x fulfills the two requirements that Xo-^i^ = “Amcti/ 
and Xofc!^ falls off fast enough at infinite spatial distances 
so that a certain surface terms vanish. One can there¬ 
fore use the Belinfante tensor T instead of the canonical 
energy-momentum tensor 0 and write 

= (A13) 

47r 

Since F is antisymmetric, the relation = 0 follows 

immediately. Noting that da{F^'*A^) = — V((^E) and 
a^(F*‘"A0) = dticpBy - (V X {(PB))^ one finds 

T™ = ^ (E^ + B2) , (A14) 

r" = (j^-<(.j)* + ^(ExB)b (A15) 


In the non-relativistic limit we may neglect A^ (so that 
E-'- —)• 0 and B -T- 0). In this case one obtains 

= ^Vip*Vip+ (A17) 

2 m Stt 

This is the expression for the energy density that will 
form the starting point for our considerations in the main 
text. 


Appendix B: Contractions rules 


We state here the contraction rules for Gaussian av¬ 
erages with the action Sq of Eq. (371. To begin with, 
the matrix P can be represented as a matrix in Keldysh 
space a^ 


0 df^,ir) \ 

dlfr) 0 ) ’ 


(Bl) 


where are two hermitian matrices. The elementary 
contraction derived from Sq reads 


{dtF,eteMKs;e3eA-^l)) 

2 

— 1^(91: ^)^q,qi ^£1 5 (-^2) 

TTVo 


where w = ei - £ 2 , ^q.qi = ( 27 r)‘^(I(q - qi), Se^,e2 = 
27r(5(ei — 62 ), and a, /3, 7 , S are spin indices. Th e dif fuson 
V was defined in Eq. (38). Starting from Eq. (B2|, one 












18 


can formulate two convenient contraction rules for the 
matrices P. The first rule can be used when the two 
matrices P appear under two different traces 


tr 


^^£ 162 (ri) tr BP^^^^{r 2 ) 


(B3) 


=-tr 

TTVq 


^'^n£i£2(i’i - r2)B^ 
Here, we denoted = i(A — asAas), and 


^ei,£4'^e2.£3' 


n£+t£-t(q) Q 

contains both the advanced and the retarded diffusions 
V and V, respectively. The following second contraction 
rule is useful when the two matrices P stand under the 
same trace 


(tr [AP^^^^{ri)BP^^^^{r 2 )]) 

— (tT[Afl^^^^{ri - r2)]tr[H] 

TTI/q \ 


(B5) 


-tr[H(T3nejeJri - r2)]tr[B(T3]^ 


'^£1,£4^£2,E3 ■ 


Appendix C: A list of contributions to Xnn" and 

In this Appendix, we provide details for the calcula¬ 
tion of Xfcfel xtini- These represent the contribu¬ 
tions of different diagrams to the dynamical parts of the 
heat-density heat-density and density-density correlation 
functions, respectively. In particular, the Appendix con¬ 
tains the definitions of the logarithmic integrals A, 
and Ji- A classification of the different types of logarith¬ 
mic integrals was discussed in Sec. |VIC[ 


1. Horizontal diagrams 


expressions stated below. A common characteristic of 
the terms corresponding to the diagrams of Fig. is that 
they contain two diffusions and, in the case of 
also the factor e^. The latter fact is the main distinction 
from the vertical diagrams to be discussed below. 

The result for the horizontal diagrams before expan¬ 
sion in Dc^ and w reads 


"Kkk,l 

Xnn,l 


dyn 

q.cj 


=27rz/o^ 


2 

q,a; 



i—a—d 


(C2) 


where 


A’m=- 


6 


X,h = - 


j i^ei -^£2 X Pe2 + i' •^£1 — 1^)^,!2J 

t/k.i/ 

2,,[D(k2+q2)-*(z. + ..)] (C3) 

t/k,i/ 

T(~T£l-|-i^ ~ ■^£2-1')^,1/] ^k+q,i/+LJ (C4) 

2 


Xic+Xid — 


- J^S2 + ^Z2+v - ^ex-v) 


/k.i/ 


X (C5) 

After expansion in and w one arrives at the simpli¬ 
fied expressions 


Xfcfe.i 

X.nn,l 


dyn 




J q,u; 


[ 


L 

1 


Ae,c.Ti(£,q,w), 

(C 6 ) 


where 

Ti =- - iuj)h{e) -I- Dc^Ioie) - iujh{e) 

+ Ti(q, w,e). (C7) 

The logarithmic integrals A are defined as 


Here, we consider contributions to the correlation 
functions originating from the expressions Xee” 

Xnn"- We specialize on those terms, whose diagram¬ 
matic representation contains a horizontal interaction 
line. They are depicted in Fig. A few remarks 
concerning these terms are in order here: diagram 

(а) contains a Hikami box. The interaction part of 
Srj=Q enters in the form (Tr[ 0 ( 73 P]Tr[ 0 cr 3 P]),)(,, diagram 

( б ) contains (Tr[(/)(j 3 P^]Tr[(/)(j 3 P^]),^, and diagram (c) 
(Tr[(/)cr 3 P]Tr[(/)(T 3 P^]),^. Here, (j) = uo (j)ou (for the sake 
of notational simplicity matrices in Keldysh spaces are 
denoted without the hat symbol here), and we used the 
notation 


(Cl) 

Each diagram displayed in Fig. I^has a symmetry-related 
partner that is not displayed, but will be included in the 


h{e) = / (P,+. - P,_.) VI^VY (C 8 ) 

D Jk,v 

lD{e) = ^ f (-?-£+. - Pe-.) D\^Vl,Vl„ (C9) 
h{e) = 11 [dePe+u + ReCif,. (CIO) 


Momenta k and frequencies v in these integrals fulfill the 
conditions < l/r and \v\ < l/r, i.e., they are con¬ 
fined to the diffusive regime. Due to the presence of A^^,^ 
in Eq. (C 6 ) important values of |£| under the integral are 
smaller or of the order of the temperature T. All three 
integrals Ii,Io,Iz are proportional to the dimensionless 
resistance p = {4:Tr'^i'QD)~^, the small parameter of the 
theory. 

Next, we turn to the logarithmic divergencies arising 
from these integrals. Important momenta for the k- 
integral in Ii lie in the range |i/|/(Pks) < k < 
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One obtains 


AW = 


f ^S + ^Z—h 


log 


DkI 


log; 


DkI 


— log -- 

6 max(|£|,r)r max(|£|,r)r 


-.(Oil) 


As will be discussed below, the double-logarithmic diver¬ 
gence arising from I\ cancels from the density-density 
and the heat-density heat-density correlation function 
after taking all corrections into account. For the inte¬ 
grals Id and Iz, relevant momenta lie in the interval 
\v\ < < I/t, and one finds 


‘A s+y ^z—u 


p\og 






1 


max(|e|,T)r’ 


(C12) 


max(|£|, T)r' 


In contrast to the three /-terms in Eq. (C7), the term 
Ji does not vanish in the limit (q, w) —)■ 0. With the 
accuracy relevant for the present calculation, it can be 
written as 


Ji(q,w,e) = 




77k-{-q,iy+a;IniV^^^. 


(CIS) 


The quantity T™ associated with the density-density 
correlation function reads 


77*" = *a;/^(£) + J 2 (q,w,e), 


(C16) 


where J 2 = —Ji and Ji was defined in Eq. (CIS). The 
quantity associated with the heat density-heat den¬ 
sity correlation function contains additional terms 

j-kk ^j-nu _ pq 2 _ + Dq^Iregie) - iujhie) 


+ 72(£). 


(C17) 


Let us comment on the appearance of these additional 
terms. If we denote the frequency associated with the 
right frequency vertex as e, then the left frequency ver¬ 
tex carries the frequency £ ± 12 due to the finite frequency 
transfer v flowing through the interaction line. Corre¬ 


spondingly, we can decompose the expression for Xi ±^2 
into a part that contains the factor £^ and a second one 
that contains the factor ev. The former part is accounted 
for by 7^"", the latter part gives rise to the additional 
terms in Eq. (C17). The integrals I 2 and I 2 are defined 
as 


hie) =f u d,iTe+. - ReVk^,, 

(C18) 

l2h) = V ^ IniRk"".- (C19) 

£ 7k,I- 


The subsequent integration in the electronic frequen¬ 
cies £ is controlled by the window function « 

ujdeJ-. Therefore, in order to obtain Ji(q, w), the in¬ 
tegral Ji(q, a;,£) may be evaluated at |£| « T. The same 
remark applies to J 2 (q, w) as well as the integrals J 2 and 
J3. We will return to the discussion of the J-terms later 
in Appendix [D| 


2. Vertical diagrams 


The integral Iregie) is not logarithmic (regular) and just 
listed for completeness. 

Let us discuss the integrals I 2 and /^ one by one. Rele¬ 
vant momenta in the integral I 2 are confined to the range 
hi < < 1 /t and one finds with logarithmic accuracy 


TT f 1 

hie) « -—p / vdeiTe+u - IFe-u) log y-p 
4e \vh 


P\q„ - \ , 

2 max(|£|, T)r 


(C20) 


Here, we discuss the contribution corresponding to the 
diagram displayed in Fig. and its symmetric partner. 
Technically, their origin is the same as for diagram 1(I>), 
i.e., the use of (Tr[(/)(j 3 P^]Tr[(/)(j 3 P^]),^ in the expressions 
for and xl^n ■ Th® bare result reads 


The integral /^ is the first in a series of (ii)-type integrals 
that will be considered. As was explained in Sec. |VI C 3| 
these integrals are mostly determined by momenta in the 
interval |z/|/(Pks) < k < ^J\v\/D, while the integration 
over the frequency v is severely limited by the combina¬ 
tion Te+v + he-v As a result one gets 


Xkk,2 

X.nn,2 


• dyn 


=2nuoVl^ 7 ^ A,,,,Pk+q..+c. 

(C14) 


TTp 


7t- t- m 1 F)k1 

ilFe+u + ^e-i^) log -y-p « P log 


max(|£|, T) 
(C21) 


Finally, let us state the integral 


An expansion up to first order in Do^ and lo gives 

Xkk,2 
Xnn,2 

(C15) 


= —2TTivQD, 


q,u; 


e‘^Tf’"ie,ci,w) 

r 2 ™(£,q,w) 


hie) [ i' ih+. + h-u)V,,.. (C22) 

e 7k, 1 / 

This integral J 2 will be discussed further in Appendix [pj 
Note that this term, unlike Ji and J 2 , does not depend 
on q and ui. 
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3. Drag diagrams 


T he tw o classes of diagrams discussed in Sec. |C 1| and 
Sec. C 2 contain a single (screened) interaction line. In 
this section we will discuss so-called drag diagrams, see 
Fig.i which form a subclass of those diagrams with two 
(screened) interaction lines. The diagrams are generated 
from 1 ^ expressions for and given in Eqs. ([33|) 
and 


It turns out that the drag diagrams do not 
contribute to the dynamical density-density correlation 
function, 


= 0 . 

Ann,3 


(C23) 


The full result for the drag contribution to the heat- 
density heat-density correlation function reads 

f f v{v - w) 

J £ J k,i/ 

X T’k,l-(25k,y -I-2?k-q,iy-u;)- (C24) 

The following two identities were used to obtain this re¬ 
sult 


[ 

e 

L 

1 


[A: 


J' en ^‘F: 


e2+i^J “ 


viv — Oj) 

0 

(C25) 


Here , £ 1, 2 = e ± w/2. The fermion frequency e in 
Eq. (C241 is associated with the right loop of the drag di¬ 


agram, while the integration in the above identities runs 
over the fermion frequency associated with the left loop. 
Upon expansion in and oj one finds 

xtT^s = j^A^^^e^Ti%e,q,uj), (C26) 


where 


r 3 “(e, q, c.) = l!^Dq^ - I^IUJ + Ue). 


(C27) 


The integrals I 2 and are rather comp licate d expres¬ 
sions resulting from the expansion of Eq. (C24) in lo and 


Dc^ and we refrain from displaying them here. With 
logarithmic accuracy, one finds 


rh _ ^ fh 

h - 2 ^ 2 , 


jh _ jh 

^2 — ^2- 


(C28) 


Relevant momenta in these integrals lie in the interval 
\v\/{Dk,s) < k < z^\u\/D, and they originate from ener¬ 
gies smaller than temperature. 

Among other terms, there also appears a finite piece 


in the expression for T 3 , Eq. (C27), 


^3(^) — 


2iyn 


lk,i> 


“(■F z + iy + -F J_^) 




(C29) 


which will be discussed in Appendix [P] together with the 
related terms Ji, J 2 , and J 2 . 


4. Regular vertex corrections 

The terms considered in this section are obtaine d fr om 
the expressions for xfr [Eq- @] and xtT [Eq- by 
taking into account nonlinear terms in the expansion of 
SQ in P-modes at the vertices. The corresponding di¬ 
agrams are displayed in Figs. [^ an d [T 0 | As it turns 
out, the drag-type diagrams of Fig. |10| vanish both for 
the density-density correlation function and for the heat- 
density heat-density correlation function. 

The expression corresponding to diagram 4a, see Fig.[^ 
reads 

xS” 4 a(q>w) = - (C30) 

X [ (-F zi-iy ~ ^ti+v — 

J k,i/ 

Due to a cancellation between the term corresponding 
to the horizontal diagram (4b) and the vertical diagram 
(4c) it is convenient to state the sum: 

xtlxbi^, w) + w) (C31) 

—7r//o^q,a; / / ['^(^ si—y £(£■ H- s2 + ix 

J £ J k,!/* 

— ^e,uj\Dk,v'Dk+<i,u+ujV^^,j. 

Correspondingly, for the density-density correlation func¬ 
tion we get the somewhat simpler expressions 

X™ 4 a(q, w) = - (C32) 

X /" (-A si-i, — T 52+1/ — 

xtT,4bi^’ ^ (C33) 

^ [ \Fe\-v — ^ei + v ~ Ae + ]T’k+25k+q,j/+<.cjf4^,/- 

Jk,ix 

Unlike for the contributions considered in the previous 
sections, no expansion in Dc^ and w is needed for the 
vertex corrections for both x!^n we may safely 

put q —)■ 0 , w —)■ 0 . 

The result can be written in the following form 


Xfcfe.4 

\.nn,A 


dyn 


= -27rwo25q.u; 


J q,a; 


e^T^’^{e,q,u) 
r 4 ™(e,q,a;) 

(C34) 


where 


rr = h{e), 
Tt = hie) + 


(C35) 

(C36) 
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The logarithmic integral /i(e) was defined in Eq. (C8) 
and the new integral is 


l 4 is) = -- [ . (C37) 

£ Jk,u 

Note that the same corrections also originate from the 
corresponding diagrams for the left vertex. Only 
is relevant and one finds the same integral as for the 
vertical diagrams: = I^- 


5. Anomalous vertex corrections 

We refer to those vertex corrections that originate from 
and Xv^ anomalous. For an illustration, see 
Fig. |11[ no anomalous vertex corrections exist for the 
density-density correlation function. The analytical ex¬ 
pressions are 


Xfefc,5a(*l ’J J ^€ 2 + 1 ^) 

X (C38) 

Xfcfc,5b(*l) f s^e,uj f {J'ei-iy -^€2 + ^) 

X (C39) 

As we are dealing with vertex corrections, we may safely 
set (q, w) —t 0. When combining these two results, one 
finds 


xS”5(q>‘^) = -27ri^o*2?q.c,; ^ l^e,u,e^Ti^{£) (C40) 

with 

Tt = -his). (C41) 

The logarithmic integral I 5 is defined as 

his) = ^l [he+u + Te-u)V^,.ReV^,. (C42) 

2 e Jk.^ 


The main contribution originates from large momenta 
D\^^ > V and one easily finds 


his) = I log 


1 

max{|£|, T}r' 


(C43) 


Another anomalous vertex correction arises from the 
drag-type diagrams of Fig. 12 (note that only the correc¬ 
tion to the 7 i vertex is written here) 


Je 

X V^ijV^q.v-ui'^k.ui'Dk,^ + T^k-q.v-Lj) 

(C44) 


Putting (q, w) —>■ 0, we find 

xtT^ei^,^^) = -27ri/oiT>q_<^ is), (C45) 

where 

rt = -I^is), (C46) 

I^ie) = + J-e-.)Fk'?.^kt2?k..ReI?k... 

S Jk,v 

(C47) 

Using the relation 

- 2i/oi^Uk^.Uk'^.ReI?k,. = (C48) 

one can transform the integral to 

4(e) = f + J-,_,)I?k.JmUk4. (C49) 

Jk.v 

Since relevant momenta are in the interval \v\/iDKs) < 
k < ./Wd, one finds with logarithmic accuracy Jg = 

1 jh 
2^2 • 

Unlike li, which are determined by very different in¬ 
tegrals, see Sec. |VIC] for the general classification of the 
logarithmic integrals, all the integrals reduce to the 
same expression. 


Appendix D: Cancellation of finite J-terms and the 
collision integral 

We reinterpret in the language of kinetics the cancel¬ 
lation of the J-terms between horizontal and vertical di¬ 
agrams for the density-density correlation function, and 
between vertical, horizontal and drag diagrams for the 
heat-density heat-density correlation function. 


1. The case of the density-density correlation 
function 

Consider the deviation of the density from its equilib¬ 
rium value caused by a weak, slowly varying external po¬ 
tential. In linear response, this deviation is characterized 
by the density-density correlation function. Our goal is 
to compare the J-terms arising during the calculation of 
the density-density correlation function to the Coulo mb 
collision integral, which is well known and read^^ZMl 

houis,x) = -2f RePk.JmUk^ (Dl) 

Jlt-M 

X [1 - Fe-uix)Feix) - S,.(2:)(J'e(a:) - Fe-^ix))] , 
where 


B.ix) 


TT 

jy 


(1 - Fe'ix)Fe'-vix)) . 


(D2) 
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This expression for Bv{x) is a generalization of the well 
known relation connecting the bosonic and fermionic 
equilibrium distribution functions, = ^/^(l ~ 

^s-iy)- 

In equilibrium, i.e., for J>(x) —>■ the collision inte¬ 
gral vanishes identically. Writing J^e(x) = J's + SJFe(x), 
the linearized collision integral reads 

Slaoids, x) =-2 [ [STe-Ax){B, - T,) 

-5Fe{x){B, + Fe-.) + 5B,{x){Fe-, - F ^)], (D3) 

where 


calculation therefore results from particle number con¬ 
servation expressed though = 0,uj) = 0. In 

the context of the diagrammatic calculation, it manifests 
itself in the absence of a mass of the diffusion. 

As a final remark on this topic let us note that the 
separation into horizontal and vertical diagrams does not 
correspond to the separation into out- and in-terms in the 
studied collision integral. 


2. The case of the heat density-heat density 
correlation fnnction 


SB,{x) = -- f SF,,{x){Fe'+,+F,,.,). (D4) 

^ Je' 

In the language of kinetics the conservation of the number 
of particles requires the vanishing of ^ SIcoii(£,x). We 
will explain here that the cancellation of J-terms origi¬ 
nating from horizontal and vertical diagrams is a result 
of the relation ^ 6Icoii{£, x) = 0. 

In an iterative approach to the kinetic problem, 
which corresponds to our perturbative treatment of the 
screened Coulomb interaction, we next use the change in 
the distribution function calculated in the absence of in¬ 
teractions, 5Fs'^\ as a zeroth order solution. It is easy to 
see that then SFe^^ is proportional to the window func- 
therefore, the bosonic distribution 


tion and that 


function remains unchanged, SB^\x) = 0. 

It is now possible to establish a connection of SIcoii 
with the diagrammatic calculation. The term propor¬ 


tional to 6Fe^^ in the linearized collision integral (D3) 
evaluated for SF = 6F^^^ is related to the horizontal dia¬ 
grams considered in Sec. |C 1[ In a similar vein, the term 
proportional to SF^^^ is related to the vertical diagrams, 

see Sec. C 2 Finally, the vanishing of 5B^\x) is directly 
related to the absence of drag-type corrections for 

Turning more specifically to the question of number 
conservation, we next focus on the expression for 
For the purpose of comparison, we reproduce here the 
integral Ji, (CI3), in the limit (q, w) —)■ 0, which arises 
from the horizontal diagrams considered in Sec. |C I| 


Ji = 


/k.i/ L 


Bu - - Fe-v) 


RePk. (D5) 


Now note that the term containing 5F^{x) in Eq. (D3) 
is proportional to Ji. In order to see this clearly, one 
should symmetrize the expression in v, using the odd¬ 
ness of A similar operation has to be 

performed for the term generated by 6F^2i, in order to 
see that it is related to J 2 = —Ji- This operation consists 
of a frequency shift e ^ e + v under the integral and 
subsequent symmetrization in ly. In summary, one finds 
that horizontal and vertical terms in the collision integral 
are determined by Ji of Eq. (D5) and J 2 = —Ji, respec¬ 
tively. The cancellation observed in the diagrammatic 


Let SFe^^ now be the perturbation caused by a 
smoothly varying gravitational potential calculated in 
the absence of interactions. In this case one finds SF^^^ oc 
£:Ae_^. As before, in terms with SFj^\, are related 


to the vertical diagrams, and terms with <5Fs^^ to the 
horizontal diagrams. Unlike for the density-density cor¬ 
relation function, however, 5B^\ which is related to the 
drag diagrams, does not vanish. 

As we will explain in the remainder of this section, 
the cancellation of constant terms between horizontal, 
vertical and drag diagrams for the heat density-heat 
density correlation function is a result of the relation 
fj. ^£SIcoii(£,x) =0. As a first step, one obtains the 


following relation with the help of Eq. (D3) and after 
shifting e ^ s + ly in the expression containing 5F^2y'- 

[ = -2 / RePk.JmUk''. 

J £ J k,i/ 

y- j[{£ + i')5F^''\B, - F,+,) - £SF^°\B, + Fe-.) 


+ eSBl°\F,_, - F,)]. 


(D6) 


The x-dependence of SF^°'> and was suppressed for 
the sake of brevity. 

We see that the terms proportional to e cancel between 
the first and the second term upon symmetrization in 
ly. This is the cancellation between the horizontal and 
the vertical diagrams encountered before for the density- 
density correlation function. The B,^ term cancels be¬ 
cause it is antisymmetric in ly. Further, we can use the 
identity siFe-v — F^) = —v^ 121: for the last term. Af¬ 
ter symmetrizing the remaining integrand with respect 
to v —v, one obtains 




/k,i^ 


RePk.JmVk^, 




(D7) 


= 0 . 


The second equality in (D7| becomes obvious upon sub¬ 
stituting the expression for 6B ^'^, compare Eq. (D4). The 


first of the two terms in the integral displayed in Eq. (D7) 
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corresponds to J 2 , i-e., it originates from the vertical di¬ 
agrams, the second term corresponds to the contribution 
from the drag diagram, i.e , to J 3 . To see this clearly, one 
should insert the identity (C48) into the definition of J 3 , 
Eq. (C29), which becomes 


This concludes our discussion of the cancellation of con¬ 
stant terms for the calculation of the heat density-heat 
density correlation function. 


J 3 = - ^ = -h- (D 8 ) 
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